Topological Scenario for Stripe Formation in Manganese Oxides 
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The spin-charge-orbital complex structures of manganites are studied using topological concepts. 
The key quantity is the "winding number" w associated with the Berry-phase connection of an Cg 
electron parallel-transported through Jahn- Teller centers, along zigzag one-dimensional paths in an 
antiferromagnetic environment of t2g spins. From these concepts, it is shown that the "bi-stripe" 
and "Wigner-crystal" states observed experimentally have different w's. Predictions for the spin 
structure of the charge-ordered states for heavily doped manganites are discussed. 
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The curious static patterns in the spin, charge, and 
orbital densities observed in manganites are currently 
attracting much attention Q. In Lai-a^CaajMnOa, the 
CE-type antiferromagnetism (AFM) has been observed 
at x=l/2 since the 1950's but recently a similar 
structure has been proposed at x=2/3 based on neutron 
diffraction experiments In these AFM structures, 
the t2g spins align in parallel along zigzag-shaped one- 
dimensional (ID) paths in the a-axis direction, while they 
are antiparallel across these paths, which are stacked in 
the 6-axis direction. This spin arrangement here is called 
the "stripe- AFM" (S-AFM) structure. 

Interestingly enough, the charge and orbital ordering 
(COO) observed experimentally is always concomitant 
with this S-AFM phase. At x^l/2, the COO has been 
confirmed by the synchrotron X-ray diffraction experi- 
ment [Q. At x—2/3, however, two different COO pat- 
terns have been reported; the "bi-stripe" (BS) structure 
in which the main building block of the COO pat- 
tern at x=l/2 persists even at x=2/3, and the "Wigner- 
crystal" (WC) structure [|, in which the COO occurs 
with the distance between charges maximized. The ap- 
pearance of the two different structures indicates that (1) 
the corresponding energies are very close to each other, 
namely, the ground state has (quasi-)degeneracy and that 
(2) the conversion between them is prohibited by a large 
energy barrier. Under these circumstances, it is of limited 
relevance the determination of which of the two states is 
better energetically based on some model Hamiltonian. 

In such a subtle situation, we focus on the origin of the 
near BS-WC degeneracy, rather than calculating which 
one is the true ground state. We follow the strategy of 
labeling these (quasi-)degenerate ground states in terms 
of a physically-motivated quantity which does not neces- 
sarily manifest itself in the Hamiltonian H. Specifically 
we focus our attention on the important role of the ID 
conducting zigzag paths in the a-b basal plane, and con- 
sidered parallel-transport of an eg electron along these 
paths through the Jahn- Teller (JT) centers composed of 



MnOg octahedra. The transport invokes the Berry-phase 
connection and we can introduce the "winding number" 
w as a direct consequence of topological invariance which 
should be conserved irrespective of the details of H. 

In this Letter, we propose that such a topological in- 
variance is a key concept to understand the complex 
states of manganites since we observe that w is always a 
good index to label the (quasi-) degenerate S-AFM states 
for x>l/2. In fact, it is found that if the energies for 
various paths considered here are plotted, its distribu- 
tion contains a multifold-band structure indexed by w. 
The observed WC and BS structures belong to two topo- 
logically different classes, characterized by the ID paths 
with w—1 and w—x/{l—x), respectively. The conversion 
between the WC and BS states is hindered by a large 
energy barrier due to this topological difference. 

Consider eg electrons coupled to both localized t2g 
spins, with the Hund's rule coupling Jh, and JT dis- 
tortions of the MnOe octahedra. Since it is the largest 
characteristic energy among those considered here, Jh is 
taken to be infinite for simplicity. This implies that the 
spin of each Cg electron at a Mn site aligns completely in 
parallel with the direction of the t2g spins at the same 
site. Thus, the spin degrees of freedom are effectively lost 
for the Cg electrons, and the spin index will be dropped 
hereafter. Since experimentally it is known that the t2g 
spins are antiparallel along the c-axis, we can assume that 
the Cg electrons can move only in the a-b plane. 

The above situation is well described in terms of H as 
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with T^i=dlidM+ dl^dai and Tzi=dl.dai- dl-du, where dai 
{dbi) is an annihilation operator for an eg electron in the 
dx2^2 (d^^2_^2) orbital at site i, a is the vector connecting 
nearest-neighbor sites, and t^^, is the hopping amplitude 
i given by i-„=-V3iJ,=- V3t?,- 3i?,-3</4 for a=x 
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and ty^=V3tl,=V3tl^= 3tl,=3t/4 for a=y §. The sec- 
ond term with the energy J' represents the AFM couphng 
between nearest-neighbor classical t2g spins normalized 
to |Si| = l. The third term describes the coupling of an eg 



electron with the {x —y )- and (3z 



-type JT modes 



(dimensionless), given by g2i and gsi, respectively. This 
term is characterized by the static JT energy i?jT Q • 

Intuitively, it can be understood that the competition 
between kinetic and magnetic energies can produce an S- 
AFM state: The system with J'=0 is a two-dimensional 
(2D) ferromagnetic (FM) metal to optimize the kinetic 
energy of eg electrons, while it becomes a 2D AFM insu- 
lator at J' to exploit the magnetic energy of the t2g 
spins. For smaller but nonzero values of J', there occurs 
a mixture of FM and AFM states one example of which 
is the CE-type AFM structure at x=l/2, schematically 
shown in Fig. 1(a). In this S-AFM state, a ID conducting 
path can be defined by connecting nearest-neighbor sites 
with parallel t2g spins. Note, however, that the shape 
of the optimal path is not obvious. A path with a large 
stabilization energy is needed to construct a stable 2D 
structure. Thus, our purpose here is to specify the shapes 
of (quasi-) stable ID paths in the S-AFM manifold. 

Let us start with the case of Ejt=0 and no electron 
correlation, which allows us to illustrate clearly the im- 
portance of topology in the present problem. In Fig. 1(b), 
the energies per site analytically obtained for various 
states are plotted as a function of J': In a 2D FM metal, 
the kinetic-energy gain is reduced by 2 J', due to the 
loss of magnetic energy per site. In a 2D AFM insula- 
tor, the magnetic-energy gain per site is 2J'. In S-AFM 
states, the optimized periodicity M for a ID path along 
the a-axis direction is numerically found to be given by 
M=2/n, in agreement with Ref. where n{=l—x) is 
the Cg-electron number per site. Thus, M is set as 4 
at x—1/2, leading to 2** paths, categorized into the four 
types shown in Fig. 1(b). The energy corresponding to 
each path is calculated and it is found that the S-AFM 
state with zigzagS path is stabilized for 0.1t^J'^0.35t. 
In the rest of the paper, this energy diagrams will not 
be shown, but a "window" for a stabilized S-AFM state 
always exists for J' of the order of O.lt at other densi- 
ties and for E'jtt^O. As shown in Fig. 1(c), the straight 
and zigzagl paths lead to a metal, while the zigzag2 and 
zigzagS paths induce a band-insulator. Since it has a 
larger bandgap (equal to t), the zigzagS path is stabi- 
lized at a;=l/2. 

The S-AFM structure with zigzagS path is nothing but 
the well-known CE-type AFM state, and our analysis 
predicts that this state is very stable at x=l/2. However, 
the same analysis for a;>2/S does not lead to a zigzag 
path as the optimized one but a straight one, which dis- 
agrees with experiment. Thus, it is necessary to find a 
quantity other than the energy to discuss the possible 
preferred paths that may arise from a full calculation 
including nonzero Ejt and Coulomb interactions. Re- 



considering the results at a;=l/2 lead us to the idea that 
the number of vertices along the path, A^v, rnay provide 
the key difference among paths. A confirmation of this 
idea is given by the calculation of energies for the 2^ and 
2^ paths at a;=2/3 (M=6) and S/4 (M=8), respectively 
As shown in Fig. 1(d), the energies can be grouped in 
(Af/2-f l)-fold bands, each of which is characterized by 
Nv{—0, 2, • • •, M). This suggests that N^, a topological 
feature, is relevant for the physics of the S-AFM states. 
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FIG. 1. (a) Schematic view of the S-AFM structure in the 
a-b plane with oxygens at the corners of squares. Solid and 
open circles represent, respectively, up and down t2g spins. 
Solid lines indicate the hopping paths for eg electrons, (b) 
Energy per site vs. J' for several magnetic arrangements at 
x=l/2 and i5jT=0. Four types of paths in the S-AFM state 
are shown, (c) Dispersion curves for the S-AFM states with 
straight, zigzagl, zigzag2, and zigzagS path. Here ao is the 
lattice constant in the a-b plane, (d) Energies per site for 
all possible paths at x=2/3 and 3/4 without JT distortions. 
Circles, squares, diamonds, up-triangles, and down-triangles 
denote the results for Nv=0, 2, 4, 6, and 8, respectively. 

Let us include the JT distortion to substantiate our 
ideas. By writing the JT modes in polar coordinates as 
921 =91 sin 6*1 and (73i=(7iCOS 6'i, "phase-dressed operators", 
dai and dbi, are introduced as dai = e'^'/^[dai cos(6'i/2) -|- 
d6isin(6li/2)] and dM = e''^'^\-daiSini0i/2) + dbiCos{0i/2)] 
with e*®'/^ representing the molecular Aharonov-Bohm 
effect. The amplitude qi is determined by a mean- 
field approximation |jl^, while the phases, ^i's, are in- 
terrelated through the Berry-phase connection to pro- 
vide the winding number w along the ID path as 
w—§(jdv-'\/6/{2Ti) ||ll[, where C forms a closed loop for 
the periodic-lattice boundary conditions fl^ . 

Mathematically w, a Chern number, is proven to be 
an integer pc| ]. In this system, it may be decomposed 
into two terms as w—Wg-\-wt. The former, Wg, is the 
geometric term, which becomes (1) corresponding to 
the periodic (anti-periodic) boundary condition in the 
Cg-electron wavefunction. The discussion on the kinetic 
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energy leads us to conclude that the state with w^=0 
has lower energy than that with Wg=l for a; > 1/2 p3[ , 
in agreement with the two-site analysis [0. Thus, Wg is 
taken as zero hereafter. 




FIG. 2. (a) A typical building block for a ID path for 
an eg electron with JT distortions, (b) General structure of 
the lowest-energy-state path and the renormalization scheme 
for the vertices a and The thick (thin) line denotes the 
straight-line part with (without) an Cg electron localized on 
it. The solid circle and diamond denote, respectively, the 
bare vertices, a and /3, while open up- and down-triangles 
indicate the renormalized vertices, a and /3. Note that the 
periodicity of the ID path is given by M=2/n=2/(l — s). (c) 
Groups of ID paths derived from mother states with m=0 
and 1. Paths in the first column corresponding to the mother 
WC structures with w—2m+l, which produce daughter states 
with u;=2m-|-2, 2m+3, ■ ■ •. 

To show that only A^v determines the topological term 
Wt, let us consider the transfer of a single Cg electron 
along the path shown in Fig. 2(a). On the straight-line 
part in the x-(?/-)direction, the phase is fixed at 0x=27r/3 
{9y—4:Tr/3), because the e^-electron orbital is polarized 
along the transfer direction. This effect may be called an 
"orbital double-exchange (DE)" in the sense that the or- 
bitals align along the axis direction to make the transfer 
of the electron smooth, similarly as the FM alignment of 
t2g spins in the usual DE mechanism. Thus, Wt does not 
change on the straight-line part of the path. However, 
when the electron passes the vertex a (/3), the phase 
changes from 6-^ to 9y {9y to 6'x), indicating that the 
electron picks up a phase change of 27r/3 (47r/3). Since 
these two vertices appear in pairs, Wti—w) is evaluated as 
Wt = (A^v/2)(27r/3-K47r/3)/(27r)=Afv/2. The phases at the 
vertices are assigned as an average of the phases sand- 
wiching those vertices, 0a=7r and Op=0, to keep Wg in- 
variant. Then, the phases are determined at all the sites 
once ^x, 9y, ^q, and 9p are known. 

Now we include the cooperative JT effect, important 
ingredient to determine COO patterns in the actual 
manganites. Although its microscopic treatment is in- 
volved, we can treat it phenomenologically as a constraint 
for macroscopic distortions pO| , energetically penalizing 
w=0 and M/2 paths. In fact, it is numerically found that 
w=l, 2, • • •, M/2~l paths constitute the lowest-energy 
band and they can be regarded as degenerate, since its 
bandwidth is about O.Olt, much smaller than the inter- 



band energy difference («0.1i). Summarizing, the coop- 
erative JT effect gives us two rules for the localization of 
Cg electrons; (i) they never localize at vertices; (ii) if an 
electron localizes at a certain site on one of the straight 
segments in the cc-direction, the other localizes on one of 
the straight segments in the y-direction. 

Applying these rules, we obtain a general structure for 
the lowest-energy path as shown in Fig. 2(b). Important 
features are the "renormalized vertices" , a and (3, abbre- 
viated notations to represent the set of straight- line parts 
that do not have eg electrons. The winding number as- 
signed to d {(3) is 1/3-I-u'q, (2/3-f W/j), where the number 
of vertices included in d {(3) is 1+2wq, (14-2^/3). Thus, 
the lowest-energy path is labeled by the non-negative in- 
tegers Wa and Wfj, leading to a total winding number 
w—l+Wa+wii. Although the topological argument does 
not determine the precise position at which an Cg elec- 
tron localizes in space, it is enough to regard a charged 
straight-line part as a "quasi-charge" . Since the quasi- 
charges align at equal distance in the WC structure, the 
corresponding path is labeled by Wa=wp=m^ with m 
a non-negative integer. By increasing wp keeping Wa 
fixed, we can produce any non-WC-structure paths with 
Wa—m and wp—m+l, m+2, ■ ■ ■ (see Fig. 2(c)). In this 
way, the WC structure with ix'=2to+1 can be consid- 
ered the "mother state" for all non-WC-structurc paths 
with w—2m+2, 2to-|-3, • • •, referred to as the "daughter 
states". The states belonging to different m's are la- 
beled by the same w, but a large energy barrier exists for 
the conversion among them, since an Cg electron must be 
moved through a vertex in such a process. Thus, the state 
characterized by w in the group with to, once formed, it 
cannot decay, even if it is not the lowest-energy state. 

Note that the topological argument works irrespective 
of the detail of H, since w is a conserved quantity. How- 
ever, it cannot single out the true ground state, since the 
quantitative discussion on the ground-state energy de- 
pends on the choice of H and approximations employed. 
In fact, either the BS or WC structure can be the ground 
state, but in view of the small energy difference, their rel- 
ative energy will likely change whenever a new ingredient 
is added to H. Especially, the Coulomb interactions will 
be important to decide the winner in the competition 
between these structures ]T5| . 

Now we analyze the charge and orbital arrangement in 
Lai-ajCa^^MnOs, in which the experimental appearance 
of the BS structure provides key information to specify 
the ID path. Since the quasi-charges exist in a contigu- 
ous way in the BS structure, its path is produced from 
the mother state with to=0 (see Fig. 2(c)). In particular, 
the COO pattern in the shortest ID path is uniquely de- 
termined as shown in Figs. 3(a)-(c). (To depict these fig- 
ures, we performed a mean-field calculation for EjT—2t, 
but the essential physics does not depend on either the 
approximation or the choice of i?jT-) At a;=l/2, the path 
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is characterized by 'w=l which is the basic mother state 
with m=0. The COO pattern shown in Fig. 3(a) leads 
to the CE-type AFM state |^. Those in the paths with 
w=2 and 3 are nothing but the BS structures experi- 
mentally observed at x=2/3 and 3/4 Note that the 
short-period zigzag part explains the peculiar feature ex- 
hibiting small oscillations in qi at less- (non-) distorted 
Mn''+ sites, as suggested in experiments 

It may be assumed that the long-range Coulomb inter- 
action V destabilizes the BS structure and transforms it 
to the WC structure but this is not the case; for the 
BS — > WC conversion with the help of V, an Cg electron 
must be on the vertex in the path with w—2 or 3 (see 
Figs. 3(b) and (c)). This is against rule (i) and thus, the 
BS structure, once formed, is stable due to the topologi- 
cal condition, even including a weak repulsion V. 

In the group of m=0, the WC structure appears only 
in the path with w=l. Thus, the WC-structure paths 
with w=l at x=2/3 and 3/4 are obtained by simple ex- 
tension of the straight-line part in the path at a;=l/2 
(see Figs. 3(d) and (e)). The detailed charge distribu- 
tion inside the quasi-charge segment is determined by a 
self-consistent calculation with the JT effect, leading to 
the WC structure. (A similar result can also be obtained 
for a weak V.) Even if the non-WC structure occurs for 
k;=1, it is unstable in the sense that it is easily converted 
to the WC structure, because no energy barrier exists for 
an Cg electron shift along the straight-line part. 



nontrivial charge and orbital arrangement will be stabi- 
lized in agreement with the experimental observation . 

Two comments are in order. (1) To understand the ob- 
served BS structure, we focussed on paths produced from 
the m—0 mother state, but we could have used a m 7^ 
mother state. Thus, we anticipate the existence of a rela- 
tion similar to Eq. (^) in a higher hierarchy (with much 
longer periodicity), leading to a devil's staircase struc- 
ture in the ID path. (2) The S-AFM structure with the 
zigzag path for w—1 explains the AFM phase observed 
in experiment at x=2/3, but we further predict a similar 
S-AFM state at x=3/4 in the WC structure. Our theory 
also predicts a peculiar spin pattern in the BS structure. 
Observation of such structures will provide a stringent 
test for the validity of our topological scenario. 

In summary, the BS and WC structures in the man- 
ganites have been classified using the winding number 
w associated with the Berry-phase connection along the 
zigzag ID path. Predictions are made for novel spin and 
orbital states in heavily doped manganites. 
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(b) x=2/3 (w=2) 



(c) x=3/4 (w=3) 




FIG. 3. (a) Patli with w^l at x=l/2 for EjT=2t. At each 
site, the orbital shape is shown with its size in proportion 
to the orbital density, (b) The BS-structure path with w=2 
at x=2/3. (c) The BS-structure path with 10=3 at a;=3/4. 
(d) The WC-structure path with w=l at x^2/3. (e) The 
WC-structure path with iu=l at x=3/4. 

The above topological analysis show that (1) the WC 
structure is made of ?i'wc=l zigzag paths and that (2) 
the BS one contains a shorter-period zigzag path with 
WBs^M/2~l^x/{l~x). Note that on the BS path, 
the less-distorted Mn'^^ sites occupy all the vertices {N^ 
equals the number of Mn'*+ ions), while the heavily dis- 
torted Mn''+ sites appear in pairs (the number of Mn''+ 
ions equal to 2). Thus, wbs is rewritten as 



WBS 



2 



Number of Mn''^"'' ions 



Number of Mn^+ ions 1 



(2) 



Since wbs is an integer, we can predict that at specific 
values of xl—WBs/i^+WBs)], such as 1/2, 2/3, 3/4, etc., 
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